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A new representation of Dirac's delta-distribution, based on the so-called q-exponentials, has been 
recently conjectured. We prove here that this conjecture is indeed valid. 

o : 

(N| ■ I. INTRODUCTION 

Tsallis and Jauregui have recently conjectured |7j a representation of the Dirac delta distribution, which they 
call S q (x), based on q— exponential functions. However, they could not prove their conjecture and used numerical 
t-H ■ experiments that suggest its validity. In this note, we provide a rigourous mathematical approach to this problem 
| and prove their conjecture by recourse to the notion of superstatistics. 

O II. q- EXPONENTIALS AND SUPERSTATISTICS 
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Statistical Mechanics' most notorious and renowned probability distribution is that deduced by Gibbs for the canonical 
ensemble |2|, usually referred to as the Boltzmann- Gibbs equilibrium distribution 



exp(-/3#i) 

Pg(i) = 7 , (1) 

Zbg 

with Ei the energy of the microstate labeled by i, /3 = 1/ksT the inverse temperature, fee Boltzmann's constant, and 
Zbg the partition function. The exponential term Fbg = exp (—/3E) is called the Boltzmann-Gibbs factor. Recently 
Beck and Cohen Q have advanced a generalization, called superstatistics, of this BG factor, assuming that the inverse 
temperature j3 is a stochastic variable. The generalized statistical factor Fqs is thus obtained as the multiplicative 
convolution 

O , 

F G s= r ^/(/3)exp(-/3£), (2) 
Jo P 

where /(/?) is the density probability of the inverse temperature. 
As stated above, (3 is the inverse temperature, but the integration variable may also be any convenient intensive 
parameter. Superstatistics, meaning "superposition of statistics" , takes into account fluctuations of such intensive 
■ - - 1 parameters. 

Beck and Cohen also show that if /(/?) is a Gamma distribution, nonextensive thermostatistics is obtained, which is 
of interest because this thermostatistics is today a very active field, with applications to several scientific disciplines 
4-6]. In working in a nonextensive framework, one has to deal with power-law distributions, which are certainly 
ubiquitous in physics (critical phenomena are just a conspicuous example Q). Indeed, it is well known that power- 
law distributions arise when maximizing Tsallis' information measure 
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H q (/) = — - [1 - J f(xYdx \ , (3) 

subject to appropriate constraints, where q ^ 1 is a real positive parameter called the nonextensivity index. More 
precisely, in the case of the canonical distribution, there is only one constraint, the energy E, i.e. (X 2 ) = E > 0, and 
the equilibrium canonical distribution writes in the case q > 1 

/,(*) = ^ (1 - (1 - q)^ 2 )^ , 
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where j3 q and Z q stand for the nonextensive counterparts of (3 and Zbg above 
In the rest of this paper, we'll assume as in Q that 1 < q < 2. 

l 



Let us choose the banch cut 
l 



-°°> 1-g 



along the negative real axis and define the q— exponential function for 



z e 



-°°'-T^ 



as 



e q (z) = (l + (l-q)z)— . (4) 
This allows us to rewrite the equilibrium distribution in the more natural way 

fq ( X ) = 7T e 9 (-^ 2 ) ■ 

It is a classical result that as q — > 1 + , Tsallis entropy reduces to Shannon entropy 

/+oo 
f(x)logf(x). (5) 
-oo 

Accordingly, the q— exponential function e q (x) converges to the usual exponential function e x . 

III. PROOF OF JAUREGUI-TSALLIS' CONJECTURE 
A. Definitions and Notations 

Recall the formula 



5{t) = 1. f e~ mt du. (6) 



We intend to provide a generalization of this relation; namely, we prove the following representation conjectured by 
Tsallis et al. assuming 1 < q < 2: 

S(t) = — [ e q (-iut) du. (7) 

for some constant c q . We begin by recalling the mathematical meaning of ©. 
Definition 1. A function (p is called rapidly decreasing if if is C°° and if for all integers k,£ 

lim x k ip {e) (x) = 0. 

x— »±oo 

Let S be the set of the rapidly decreasing functions on R and by S' the set of the continuous linear functional 
over S. For (p £ S, its Fourier transform J-((p) is denoted by (p. 

We know from Rudin [j| p. 184 theorem 7.4] the following 

Proposition 1. The Fourier transform J- is a continuous linear mapping of S into S. 

Definition 2. Let now f be a bounded measurable function flJl] - We let Tf be the linear continuous mapping: 



VcpeS (T f ,<p)= f(tMt)dt 



B. Proofs 

In order to prove the usual representation ((6]), we simply have to show that for all ip G S 

du (T e -,ut,(p) — 2ir<p(Q). 
Of course (T e -»««, ip) = ip(u). Hence the result. 

We now turn to the proof of ((TJ). In this respect, let us pick a <p £ S. We have 

{T eq {-iut),v) = j e q (-iut)tp(t)dt 

= J E w e~ lut{q - 1)w 'tp(t)dt 

where we have used the equality 

e q (-tut) = E w e~ mt{q ~ 1)w . (8) 

Here 

A / / x -id... -1. 



= — -r- / g(w)e- w w — - L dw 

is the expectation of g(W), where W is a Gamma distributed random variable with shape parameter -^j and g some 
function such that the above definition makes sense. 

We note that |8]) expresses the fundamental principle of the superstatistical theory. 
On the other hand, we have 
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,J w^~ 1 \cp(t)\dtdw < / \ip{t)\dt 



As obviously ip is summable, we can apply the Fubini-Lebesgue theorem and we obtain 

(T eqi - m t)^) = E W J e- mt ^ w <p(t)dt 
= E w <p(u(q-l)W) 

= — ; — : — - / e~ w w ~^ rT ~ 1 (p(u (q — 1) w)dw 

Now, consider 

/ (T eq (-zut),¥>)du. 
JR 

As q < 2, we have, by the change of variable u t— > v = u (q — 1) w, 

dw J \e- w w*^- l y(u (q - 1) w)\du = J \^)\dv J e - v , V) .^ I -2 dw < ^ 

since, by Proposition [TJ (p G S. Thanks to the Fubini-Lebesgue theorem, we deduce that 

it \j 1 i ' V{v)dv f _ w i _ 2 

r(^r - 1) 

9-1 rr2^(0) 



(9-i)r(^) 

We have proved the result with 
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IV. REMARKS ON JAUREGUI-TSALLIS' APPROACH 



In their approach [7|, the authors chose an empirical approach starting from the usual q = 1 case: expressing the 
Dirac delta as the limit 

5{x) = - lim [ +L e -^dk=- lim 

C\ L-s-+oo J _ L Cl L->+oo X 

the normalization constant c\ can be obtained formally as 

f + °° sin fix) , 
ci = 2 lim / ax = 27r. 



L— >+oo 

The extension to q— exponentials reads 



2 sin ( arctan((<7 — 1) TV 

5 g (x) = — r — lim 



(2 - q) c g L^+oo + l)L 2 x 2 )^> 

so that the normalization constant c q can be obtained formally as 

• (§E? arctan ((<?-!)£)) 



+oo sin I 

lim / — , „ ' dx. 



2-qL^+ocJ_ OQ x + l)L 2 x 2 )^y 
This integral can be equivalently expressed, using the change of variable z = tan 9 as 



sin ^0 (cosfl)^ 1 
I g 4 2 / " ^ | „ d0. 



sin 



The authors then evaluate the integral I q for a finite number of rational values of the parameter q <= ]1, 2[ only, for 
which the symbolic computation software Maple gives the value I q = j. However, this approach can be circumvented 
given the fact that the integral I q can be found in [T3, 3.638.3], with the value I q = | Vq G ]1, 2[ so that empirically, 

C 9 — 2-g' 



V. CONCLUSION 



We have proved that the representation of the Dirac delta distribution ([7]) using q— exponential functions, as conjec- 
tured by Tsallis et al., is valid. In particular, (i) we compute the exact normalization constant in the representation 
of the Dirac delta and (ii) we explicit the set of functions for which this distribution acts as the Dirac delta. 
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